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We present a model for the description of the evolution of contacts among individuals in a network. 
At each time step each individual is associated with a domain or neighborhood of fully connected 
agents. The dynamics of this changing neighborhood will later be translated into a situation where 
the links between individuals are also dynamic. A characterization in terms of the parameters that 
govern the evolution of the network and a comparison to previous work on Small World networks is 
presented as well. 
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INTRODUCTION 



The study of networks and their appHcations within 
social sciences has in recent years been a rich source of 
interdisciplinary research. In particular, development of 
mathematical models based on simple social interactions 
formulated with graphs or networks, with certain topo- 
logical properties, has accumulated much interest. The 
concept of the small world was first introduced by Mil- 
gram Q in an attempt to describe the intrinsic prop- 
erties of social communities and the relationships that 
exist amongst the members thereof. A more mathemati- 
cal approach was later introduced to understand the un- 
derlying aspects of Milgram's networks Q. Small World 
Networks (SWN) are in essence, those networks which 
are comprised of a specified amount of random and reg- 
ular lattice connections. Driven by strong evidence that 
certain complex networks share a close resemblance with 
actual economic, social, and biological networks, a new 
field has been developing within physics Such real 
world applications have naturally led to the proposition 
of temporally evolving networks [1,1^. We will refer to 
dynamic networks as those which rewire links as time 
progresses. Networks in which the links remain static in 
time once the desired architecture has been achieved, are 
known as quenched networks. Whereas networks with 
links that change randomly at every time step are called 
annealed networks. An example of the latter, which is 
akin to the model developed here, may be found in 
The model we present consists of subnetworks which are 
only linked to each other through time. A dynamic do- 
main network (DDN) is thus a type of annealed network, 
though at a given moment the subnetworks or domains 
are not connected to each other. As explained in the 
next section, the model we propose may be associated 
with a phenomena where the only mechanism of trans- 
mission or propagation throughout the network is a type 
of diffusion. Though the model can be generalized to in- 
clude other features such as long range transmission (as 
found in SWNs), we are for now interested in analyzing 
properties due solely to this specialized diffusion. The 



possibility of analyzing different dynamic degrees of free- 
dom independently will lead to a deeper understanding 
of diffusion within networks. In the following sections 
we will consider in more detail the features guiding our 
analysis, 

• Two measurable parameters, uj and r, which may 
be associated with the SWN parameters, clusteri- 
zation and average path length. 

• Three types of domain movement, pd, Pf, and Ps, 
which are called the dynamic processes of the sys- 
tem, are analyzed 

• Numerical analysis and limited analytic calculation 
of oj and t as functions of the dynamic processes 
demonstrate common characteristics with SWNs. 

• Diffusion coefficients associated with each dynamic 
process are numerically and analytically calculated. 

• The DDN model reproduces features found in pre- 
vious work on disease propagation and further 
shows the affect of each dynamic process on the 
total population infection. 

THE MODEL 

The model is based on a two dimensional array of indi- 
viduals, each one situated on the node of a matrix of size 
n X n. The system is divided into N domains, or fully 
connected subnetworks, of size fc^. In what we call the 
regular case, all domains have the same shape and size. 
We consider that an individual is linked to the other in- 
dividuals belonging to the same domain, but not to the 
rest of the system. The boundaries that demarcate the 
domains change in time and thus may be considered as in- 
dependent neighborhoods exchanging members through 
dynamic links in time. The architecture of the underlying 
dynamic network is defined by the history of the changes 
of the boundaries. A set of three processes govern the 
behavior of these domains and thus three different pa- 
rameters will be considered. At each time step the whole 



2 



set of domains may be displaced or shifted from their 
former positions with probability pd- It is also possible 
that domains change individual size with probability Ps 
and change form or shape with probability p/. Each of 
these processes can act separately, but when all three take 
place at the same time, a superposition of their effects oc- 
cur. To characterize some of the dynamic properties of 
the evolving network we have chosen a few observable 
quantities as suggested in We consider a measur- 
able global variable analogous to the average path length 
in SWNs. Over all domains and throughout time, we 
analyze the time necessary for an 'active' state to prop- 
agate across the entire system. Initially, one individual 
is chosen to be active. Any inactive individual in the 
same domain is immediately activated by the presence 
of at least one active individual in the same neighbor- 
hood. Once activated, the individual remains as such for 
all time. We measure the time involved r^, to activate the 
node i within the system. The mean persistence or acti- 
vation time T, is calculated as defined in |i2i]. This value 
is studied for different sets of Pd, Ps, and p/. As a second 
characterization of the system, we analyze the changes 
in the composition of each individual's domain from one 
time step to another. Again, in accordance with 0, this 
local measurable variable may be associated with a clus- 
tering coefficient ■ As a general definition we call uji {t) 
the overlap or intersection of domains during m consec- 
utive time steps, where in the present work m — 2. The 
mean neighborhood overlap is then calculated for varying 
values of Pd, Ps, and pf. Compared to previous models 
of dynamic networks, this model may be associated with 
diffusive processes. This supposition will be reinforced 
with the calculation of associated diffusion constants. 



ANALYTIC CONSIDERATIONS 

As a measure of the common area of two domains be- 
tween consecutive time steps, the overlap may be written 
analytically in a general form. In the model we have con- 
sidered that the allowable sizes take on only quadratic 
values so that the square shape is always accessible, par- 
ticularly for the case when the form does not change at 
all. We thus compare all possible forms of the permit- 
ted sizes between a before domain and an after domain. 
The expression for the mean overlap for any set of values 
P = {Pd,Ps,Pf) can be written in general as a binomial 
distribution 



i=0 j=0 k=a 

a.^kima-Pdr-Ysii-Psr-'pfii-Pf) 



(1) 



The coefficients aijk{p) must be calculated by consider- 
ing the superposition for any succession of two allowed 
states of the system for a given set of values p. Since we 



preserve the rectangular shape of the domains, the over- 
lap is determined by the divisors of domain sizes in the 
before and after configurations. This allows us to write 
a counting scheme for the overlap when pd — 0. 
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(2) 

With this, the overlap is calculated for the four extreme 
situations when ps and p/ are either zero or one. The 
parameters Nm and Nk indicate the number of divisors 
of the sizes and fc^, for each before and after domain, 
respectively. The value is the maximum size allowed. 
To calculate the case when ps = Pf = 1, one must set 
a = 7m = 7fc — 1; this sums over all forms and all 
sizes bounded by n. For zero change in form and size 
Ps — Pf — 0, the parameters must be set to a = n, 
7m — ^"2""'" + 1, and "fk — + 1- Changes in size 

Ps — 1, and no change in form P^ — correspond to 
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1, and 7fc = 



1. Changes in 



form Pf = 1, and no change in size Ps — 0, correspond 
to a = n and jm = 7fe = 1- These extreme values define 
four of the eight coefficients in Eq. 

The underlying diffusive process occurring in the net- 
work as time proceeds is another interesting feature we 
would like to address. We begin by considering one par- 
ticle located within an initial domain. This test subject 
then performs a walk to neighboring domains as the do- 
main walls change and an overlap between the original 
domain and a neighboring one permit movement of the 
individual. Indeed, the process may be considered as a 
random walk, characterized by a probability distribution 
particular to each dynamic process. A brief analysis of 
what movement might be expected can be done by con- 
sidering the possibilities of change associated with each 
process. As an example, consider the case of a dynamic 
domain originally comprised of nine sites. A probability 
can be calculated for a given random walker located at 
the center of the original domain to step to the center 
of one of the neighboring domains with which it over- 
laps in a given time step. Two different jumps should be 
considered; those to neighboring domains that share a 
side with the original one, and those to any of the other 
four domains along the diagonals. We will call a the 
probability of the first type of jump, and /3 the proba- 
bility corresponding to jumps of the second type. The a 
and /? probabilities for each of the three processes can be 
calculated by considering the fractions of all the overlap 
situations that will lead to a given jump, where we have 
calculated ad = 0.12, (3d = 0.04, as = 0.024, (3s = 1/256, 
af — 1/9, (3f = 0. These values are used to evaluate 
a random walk process and thus an associated diffusion 
coefficient. 

A random walker in one dimension with equal likeliness 
to jump to the right or left is analyzed. This probability 
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to move is given by p(a + 2/5), where a and /3 corre- 
spond to the values mentioned previously and p is the 
value characterizing the probability for a change in the 
position, size, or shape. Further, there is a probability 
to stay in the same site, 7 = 1 — p(2a + 4/3). This leads 
to a calculation of the mean square displacement (a;^), of 
a set of particles equal to 18p(a + 2f3)t By consider- 
ing the appropriate values for each dynamic process, we 
find that the diffusion coefficients for each are given by 
Dd = l.Spd, D, = 0.28p, and Z?/ = pf. 

NUMERICAL RESULTS 

Extensive numerical simulations of the described 
model for networks of = 10^ to 10^ nodes and k = 3 
to 10 have been performed. A typical realization starts 
with the random election of an active individual. In the 
successive time steps, the active condition propagates 
throughout the system until the system is entirely ac- 
tivated. We calculate the mean activation time of the 
system as 

1 

i=l 

where is the activation time of the i node. We also mea- 
sure the mean neighborhood overlap between two consec- 
utive time steps as 

t i 

where the sum over time is performed over all the time 
steps until total activation at time ta- In FigQlwe plot 
the values of r and oj for the specific cases when two of 
the three p values are zero and the remaining one varies. 

A comparison between the analytic calculations of lu 
and the numerical results are displayed in Fig. |3 Again, 
only the limiting cases are considered. The lower curves 
correspond to changing either Ps or p/, while pd and the 
remaining p value are both zero. Correspondingly, the 
upper curves depict a varying ps or pf, while pd is again 
equal to zero, but the remaining term is now constantly 
one. The continuous line is obtained by analytic calcu- 
lation on top of which we have plotted the numerical 
results. 

An interesting aspect that leads us to consider diffu- 
sion within the system, is the evolution of the size of the 
cluster of active individuals. For this we have performed 
several calculations considering an initial active agent lo- 
cated at the center of the system. For several values of p 
and considering the variance of only one p parameter at 
a time, we have evaluated the mean radius of the growing 
nucleus. The radius evolves linearly with time and thus 
the velocity of propagation of the active state v, is easily 
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FIG. 1: The bold line r, and the dotted line ui as functions 
of p, where p — (squares), p = Ps (triangles), and p = Pf 
(circles). Non- varying p values are set equal to zero, r is 
normalized to the total activation time when p — 10~2. 
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FIG. 2: The average overlap lu, as a function of p, where 
p — Ps (triangles), p = Pf (circles), and Pd = in all cases 
plotted here. Solid symbols indicate that the non-varying p 
parameter is one, whereas open symbols indicate the non- 
varying p parameter is zero. 

calculated. In Fig. El f is depicted as a function of p for 
each dynamic process. The inset shows w as a function of 
the mean overlap. Apparent from the figure, the velocity 
adopts a different behavior for each of the three dynamic 
processes. 

In conjunction with analytic calculations of the diffu- 
sive process, we have performed numerical simulations 
to measure the mean square displacement (x^}, where 
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FIG. 3: The average velocity of propagation v, as a function 
of p, where p — pd (squares), p = Ps (triangles), and p = Pf 
(circles). In the inset we plot f as a function of u> 

{ ) indicates the mean value over individuals and x is 
the position of a random walker on the a; axis. From 
the analysis of this quantity throughout time, we obtain 
a value for the diffusion coefficient. Fig. 31 shows this 
value as a function of each of the three processes defined 
by Ps, Pd, and pf. To compare numerical results with 
the a and /? values obtained before, we have considered 
several individuals performing a random walk in the sys- 
tem. Walkers move along the centers of the domains. 
A given jump is allowed only if superposition between 
two domains occurs in two successive time steps. The 
results correspond to networks undergoing the different 
processes governed by the aforementioned p parameters. 
The fitting by minimum squares methods gives us the 
following values: Dd = l.Gpd, Dg = 0.27ps and Df = Pf, 
which are in agreement with previously calculated values. 

DISEASE PROPAGATION 

The DDN model that we have proposed compares to 
previously obtained results of disease propagation within 
Dynamic Small World Networks 9] in the following way. 
We consider a standard model dii] for an infectious dis- 
ease with three stages: susceptible (S), infectious (I), and 
refractory (R) . Any susceptible individual can become in- 
fected with a given probability by an infected individual 
within the same domain or neighborhood. The infec- 
tion cycle ends when the element reaches the refractory 
state after r time steps. The refractory state is a perma- 
nent condition and thus the individual cannot be infected 
again. The algorithm goes as follows. At each step an in- 
fected element i is chosen at random. If the time elapsed 
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FIG. 4: Diffusion coefficients as functions of p, where p = 
Pd (squares), p = ps (triangles), and p ~ Pf (circles). Solid 
lines indicate the analytic calculation over which the numeri- 
cal results are plotted. 



from the moment ti when it entered the infection cycle 
up to the current time t is larger than the infection time 
r, the element i becomes refractory. Otherwise, one of its 
neighbors j of i is randomly selected. If j is in the suscep- 
tible state, contagion occurs. Element j becomes infected 
and its infection time tj = t is recorded. If on the other 
hand, j is already infected or refractory, it preserves its 
state. Since each time step corresponds to the choice of 
an infected individual, the update of the time variable 
depends on the number Ni{t) of infected individuals at 
each step, t ^ t + \/Ni(t). Once there is no change 
in the number of infected individuals for the duration of 
time necessary for all infected individuals to transform 
into the refractory state r, then activation time has been 
achieved and the process stops. Here we set r = 3 for 
the infection time which insures that for intermediate val- 
ues of the p parameters, the disease spreads over a finite 
fraction of the population. Moreover, we consider an ini- 
tial condition where there is only one infected element; 
all the other elements being susceptible. The considered 
initial condition thus represents an initially localized dis- 
ease. The infection will remain localized during the ini- 
tial stages and will propagate according to the behavior 
of the domains. In Fig. 0we plot the proportion of total 
infected individuals if as it varies with different values 
of p. Two of the p values are maintained fixed while the 
remaining one varies between zero and one. 

In comparison with the results obtained in 0, we find 
that the threshold values of p for disease propagation are 
of the same order. The sharp transition from the non- 
propagative regime to the propagative one as a function 
of the structure of the network is one of the most in- 
teresting aspects presented by epidemiological models on 
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FIG. 5: The fraction of infected individuals with respect 
to the total population if, as a function of p, with p = 
Pd (squares), p = ps (triangles), and p = Pf (circles) 

networks. The maximum proportion of infected individ- 
uals shows that a fraction of individuals remain not in- 
fected. Similar values of final infected individuals have 
been obtained in 0. The non-monotonic behavior of this 
quantity has also been observed in other emdemiological 
models based on Small World Networks A possible 
explanation for the decrease of infection for large values 
of Ps involves the fact that in the DDN system, there is 
a bias towards decreasing the size of a domain. A high 
level of dynamics will temporally isolate infected individ- 
uals. It is because of this limitation of the model that 
the infection cannot propagate as efficiently for large Ps 
and a smaller number of infected individuals is detected. 

CONCLUSIONS 

If we aim to compare the system presented here with 
previous work on dynamic Small Worlds, we should limit 
this comparison to cases when shortcuts only to relatively 
close neighbors are considered. This restriction is due to 
the fact that in the present model, new links can only 
be established with nodes occupying neighbor domains. 
An interesting aspect despite these Hmitations is that we 
still observe some of the features found in dynamic SWN 
models without such restrictions. In the first part we 
have studied the behavior of the mean persistence time 
r, which is similar to the SWN path length, and mean 
neighborhood overlap uj, analogous to the SWN clusteri- 
zation. Though no analysis of how the domain size scales 



with these quantities was done, the qualitative similari- 
ties ought to be noted. That is, there exists a region in 
the p parameter space when the system displays a mix- 
ture of the two Hmiting cases, relatively high overlap and 
relatively low activation time. The value of to can be an- 
alytically obtained for some limiting cases. Further, we 
have studied the underlying diffusive process occurring 
on the network. A diffusion coefficient was numerically 
as well as analytically calculated for the three dynamic 
processes, Ps, Pf, and pd and demonstrated the diffusive 
contribution of each. As a final application and point 
of comparison with previous work, we analyze the prop- 
agation of an infection. For small p values, an initial 
infection can not propagate. At a given p value, propa- 
gation is possible and grows rapidly to a maximum in- 
fected fraction of the entire population as p increases. 
The threshold values are similar to those previously ob- 
served for a dynamic SWN. In this case, we show that 
despite the fact that there is a transition around the same 
value of p, the behavior of the system strongly depends 
on which dynamic domain process governs the dynamics 
of the system. 
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